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ABSTRACT 

Two power series inverses of the universal form 
of Kepler's equation are derived with simple recursion 
formulas for the coefficients. The first series is in 
powers of time, the second in powers of a = 1 a, where 
a is the semimajor axis. The second series is par- 
ticulary useful for nearly parabolic orbits, since 
Kepler's equation is easily solved when a = 0, i.e., 
when the motion is parabolic. 
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LIST OF SYMBOLS 


a semi-major axis of orbit 

e eccentricity of orbit 

P semilatus rectum of orbit 

r orbital position vector at time t 
r magnitude of the orbital position vector at time t 
7 o orbital position vector at time t o 
r o magnitude of the orbital position vector at time t Q 
r time derivative of r at time t 

O O 

t time 

t reference time 

O 

6 true anomaly of orbit 

ii gravitational constant times the mass of the attracting body 
v orbital velocity vector at time t 

v q orbital velocity vector at time t o 

v o magnitude of v q 
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POWER SERIES INVERSES OF 
THE UNIVERSAL FORM OF 
KEPLER'S EQUATION 


INTRODUCTION 

The universal form of Kepler's Equation is Cl] 


Acr + By + vji - fx A ( t - t Q ), 


( 1 ) 


A = 1 - ar , 

o 


B = 


V 2 

O 

M . 


a= — 

a r 


cr = /3 3 S(cl/3 2 ), 


( 2 ) 


7 = /S 2 C( a/3 2 ) . 


The S and C functions are defined by the series 


00 


soo=£ 


i =0 


(-X ) 1 

( 2i + 3)! ’ 


( 3 ) 


( 4 ) 


1 


c <»= E 

; = A 


Given the time t and the position and velocity at time t 0 , (3 is found from 
equation (1). The position and velocity can then be found at time t from 


r = ( 1 - 7/r o ) r o + ( t - t o - ) v q . 


( 6 ) 



r r 

o 


(acr-yS) r Q + (1 -y/r) 


(7) 


r = Ay - B( cur - /3) + r Q . (8) 

Chebyshev approximations of the functions S(x) and C(x) have been obtained 
[2] which reduce significantly the computation times required for their eval- 
uation when compared to that required by the series (4) and (5) . 

We have obtained a power series solution of equation (1) for (3 as a function 
oft - 1 . The procedure simultaneously produces series for a and y, which are 
needed in the updating equations (6), (7), and (8). The recursion formulas for 
the coefficients in these series are extremely simple, making possible a com- 
puter program which is not only fast but requires very few computer storage 
locations. 

If we let a = 0 in equations (2) and (3) for a and y , equation (1) becomes a 
cubic in /3: 


~ A/3 3 +\b/3 2 + r 0 /3-/u* (t - t o ) = 0. 

D 1 


(9) 


Our second series is a perturbative power series solution of (1) for (3 as a 
function of a with t fixed, i.e., 

00 

/3= Z] ^ ai> 

i = 0 


2 


where /? 0 is the solution of (9) . 

A sufficient condition for (9) to have only one real root is that the quadratic 
obtained by differentiating (9) have no real roots. 

The discriminant of this quadratic is 

D □ B 2 - 2 r A 

O 



r 2 r 2 

O O 





where p is the semilatus rectum of the orbit. Substituting 


° 1 + e cos 6 


into the expression for D gives 

ep (e + 2 cos 0 + e 2 cos 2 0) 

^ ( 1 + e co s 0) 2 

Thus the quadratic will have no real roots provided 

e + e 2 cos 2 0 + 2 cos 0^0, (i0) 

and this is a sufficient condition for (9) to have only one real root. 

Obviously (10) holds for all hyperbolic and parabolic orbits and for all 
elliptic orbits such that 1 0[ ^ 90°. For 0 > 90° (9) has three real roots for some 
elliptic orbits. We will not pursue the matter further, however, since the uni- 
versal form of Kepler’s equation is not needed for 0 > 90°. 

In reference [3] the initial position vector r o is chosen at periapsis, and 
after finding fi from the cubic equation, a correction is found by the Newton- 
Raphson procedure. 
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In reference [ 4 ] recursion formulas are derived for the coefficients of three 
power series inverses of the elliptic form of Kepler’s equation: (1) a series in 
powers of time, (2) a series in powers of eccentricity, (3) a bivariate series 
in powers of two small parameters. 


POWER SERIES IN TIME 

We define the time from epoch as r = t - t , so that from (1) we have 

Act + By + r o {3 ~ y/Hr. (11) 

Then from the identities 

^T = TT [C(x) -3S(x>] , 
dx 2x L J 


^ [l - xS(x) - 2C( x )"] , 
dx 2x L J 

which are easily established from the definitions (4) and (5), we obtain the 
equations 

* = 7 / 3 , ( 12 ) 


y - (fi ~ aa ) b 


where the dot indicates differentiation with respect to time. 


Let 


/3= y \ A ri> 

: ss n 


= 2Z -i" 1 ' 


(13) 


4 


Then 


7 = 


i =0 


CO 

4 = J2 (i + 1) yS i + 1 T 1 , 

i = 0 


cr _ 


I 

; = n 


( i + 1 ) o'i + i ri - 


00 

7 = ( i + 1 ) 7j + j r 1 . 

i = 0 

Substituting the last six equations into (11), (12) and (13) and equating co 
efficients of like powers of T , we obtain 


A cr. 


+ 1 


+ B A + 1 + r , A 


+ 1 


= o 


(i + O a i + i = ( j + 1} ^i + i A-j - 

j = 0 


1 1 

(i + 1) 7 i + 1 = " a L (j + + l CT i-j 

j = o j = o 


all valid for i > 1. 
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However , noting that 


A„ 


- a 


y 

o 


(14) 


Aj - ri - ^1 - ° - 


(15) 


A 2 



^2 


1 

“J M . cr 2 “ 0 > 


we have the following recursion formulas for i > 2, 


(16) 


S. 

j 




i " 1 

j = i 


(17) 


( i + 1 ) 7i + i = 



S. 

J 




8 . a. _ . 

j i J 


+ l 


r o Ai + ! = - ACT i + X - B ^i + 1 


If we wish merely to calculate updated values of A, cr, or y, we let 


A x = . 


6 


4 = 


BflT 2 

2 r 

O 


in equations (15) and (16). The recursion formulas (28) then produce coefficients 
such that the series become 


fi = 



Pi 


00 



i = 0 


7 = 


7 

; = n 


y. • 


POWER SERIES IN a 

Equations (2) and (3) are equivalent to [ 1 ] 

o? /2 <y - a} /2 /3 ~ sin(a 1/4 /3) t 


for a > 0 and 


( -a) 3/2 cr = s i nh ( ( ~af/3 ) - ( -a)*/3 

for a < 0 and 
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for all a . 
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a - 


UJ 

L - 


i = 0 


uu 

y = 2Z 7i ai ’ 


i = 0 


UJ 

ZZ s > • 


i = 0 


00 


= £ p, ■ 


i = 0 


Substituting these series into (18), (19), and (22) and equating coefficients of 
like powers of a, we have, for i > 1, 


i - 1 


— / 5 2 (i+— )/3. + — /3 y . - ( i + — ) a. — ^ ' p.y . _ , 

2 ° 2 1 2 ° 1 2 1 Z—i 11 1 

i = i 


(23) 


A 


-?i = ■ £ Pj s i-j j/^ i + 1) ’ 


(24) 


r B. + By . + Act. - 0 

O ~1 7 1 1 


(25) 


Multiplying (24) by -j /3 Q and adding to (23) gives cr. in terms of /3. . Sub- 
stituting that result with (24) into (25) yields a recursion formula for the /?.. 
We also substitute the series into (18) and (19). Our final result is 


p = — B , S = P> , y = pS,cr :-p 7 , 

*0 2 ° 7 ° *0 O o ^ 
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where B is the real root of 

O 

La^ 3 + -B/3 2 + r 0 - C = 0 , 
6 2 


S = 


r + 

O 


B/3 o + Ay o 


b = — p O' 
2 ° o 


f\ =b(B + |p o AyS , 


Pl = | /3i . y i = AA “ b ' = y A " 


and for i > 1, 



(i + 1) 


G = 




Vi-i 



/3. = (BF + AG) /S 


_F 
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CT i = y Q Pi ~ G • 
Pi =(i + f)£, 


We note that the divisor S in the equation for /3. can be written 

S = - A/3 2 + B/3 + r , 

2 ' o ' o o 

and this cannot be zero if we conform to the restrictions suggested in the Intro- 
duction. 

Again, if we wish to calculate updated values of /3, a, or y, we let 


F = 


ap a. , 

*0 l l 


i - 1 

- L 

j = 1 


P; - 



(1 + 1 ) 


S. = yS. - acr. 

i ' i i~i 

so that the recursion formulas produce coefficients such that the series become 


z> 

: = n 


00 



i =0 


li 



7 : 


- L 

: = a 
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